A thermodynamically consistent quasi-particle model without 
temperature-dependent infinity of the vacuum zero point energy 
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In this paper, an improved quasi-particle model is presented. Unlike the previous approach of 
establishing quasi-particle model, we introduce a classical background field (it is allowed to depend 
on the temperature) to deal with the infinity of thermal vacuum energy which exists in previous 
quasi-particle models. After taking into account the effect of this classical background field, the 
partition function of quasi-particle system can be made well-defined. Based on this and following 
the standard ensemble theory, we construct a thermodynamically consistent quasi-particle model 
without the need of any reformulation of statistical mechanics or thermodynamical consistency 
relation. As an application of our model, we employ it to the case of (2-(-l) fiavor QGP at zero 
chemical potential and finite temperature and obtain a good fit to the recent lattice simulation 
results of S. Borsanyi et al. A comparison of the result of our model with early calculations using 
other models is also presented. It is shown that our method is general and can be generalized to 
the case where the effective mass depends not only on the temperature but also on the chemical 
potential. 
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I. INTRODUCTION 



It is generally accepted that the theory describ- 
ing strongly interacting matter is Quantum Chro- 
modynamics (QCD). In QCD with two massless 
quarks, a spontaneously broken chiral symmetry 
is restored at finite temperature and zero chemical 
potential. It can be argued [l|, Q that this phase 
transition is likely second order. If this transi- 
tion is indeed second order, QCD with two quarks 
of nonzero mass has only a smooth crossover as 
a function of T. Although not yet firmly estab- 
lished, this picture is consistent with present lat- 
tice simulations and many models. Meanwhile, at 
finite temperature and zero chemical potential, the 
transition from the confined hadronic matter to 
the dcconfincd quark-gluon plasma (QGP) is un- 
derstood now as a crossover for QCD in the real 
world 0^. Because QGP describes the relevant 
features of nature under extreme conditions, for 
example, in neutron stars (at large density [6|), or 
in the early universe (at essentially vanishing net 
baryon density [7| ) , the study of thermodynamical 
properties of QGP has attracted considerable at- 
tention over the years from the theory community. 
In principle, lattice QCD provides a straightfor- 
ward way to calculate properties of QGP, and in 
particular, its equation of state (EOS). The tech- 
niques of these calculations have been developed 
greatly in the last years, even for the case of fi- 
nite chemical potentials |8l-[l3|. Among these an- 
alytical approaches, the weak-coupling expansion 
of thermodynamical quantities shows a extremely 



poor convergence for any temperature of practical 
interest [l3J. To overcome this poor convergence, 
researchers have invented many methods with a 
rigorous link to QCD, notably the resummed HLT 
scheme [ij, Il3|j ^-functional approach [1^ Il7| . 
On the other hand, many different phcnomeno- 
logical models were adopted to describe the non- 
perturbative behavior of QGP seen in lattice sim- 
ulations of QCD. Among those models, the quasi- 
particle quark-gluon plasma model with a few fit- 
ting parameters, has been widely used to repro- 
duce the properties of the QCD plasma |18[ - [28| . 
In this model, at finite temperature, instead of 
real quarks and gluons with QCD interactions, 
one can consider the system to be made up of 
non-interacting quasi-particles, i.e., quasi-quarks 
and quasi-gluons, with temperature-dependent ef- 
fective mass. The quasi-particle model was first 
proposed by Goloviznin and Satz [T^l , and later by 
Peshier et al [19J • After some time, Gorenstein and 
Yang found that this model was thermodynami- 
cally inconsistent, and remedied this flaw by re- 
formulating statistical mechanics [20|. After that, 
Bannur also proposed a new quasi-particle model 
using standard statistical mechanics and avoided 
the thermodynamical inconsistency from the en- 
ergy density rather than the pressure [27| . Further- 
more, in Ref. [291 Gardim and Steffcns showed that 
the two models proposed by Peshier and Bannur 
were two extreme limits of a general formulation. 

Over the past few years, more considerable 
progress has been made in the quasi-particle 
model. However, just as will be shown below. 



in all previous works on quasi-particlc model the 
problem of the temperature-dependent infinity of 
the vacuum zero point energy and its influences 
have not been seriously considered. In a consis- 
tent quasi-particle model, how to eliminate this 
temperature-dependent infinity of the vacuum zero 
point energy is very important. In the present 
work we try to answer this question. 

The outline of this paper is as follows. In section 
II, we show the problem existed in previous quasi- 
particle model. In sections III and IV, adopting 
the series expansion method inspired by Walecka 
to deal with the temperature-dependent infinity of 
the vacuum zero point energy, we construct a new 
thermodynamically consistent framework of quasi- 
particle model for QGP without the need of any 
reformulation of statistical mechanics or thermo- 
dynamical consistency relation. As an application 
of our model, we employ it to the case of (2-^1) 
flavor QGP at zero chemical potential and finite 
temperature to fit the recent lattice simulation re- 
sults and compare our results with those of earlier 
models. In section V, we conclude our work with 
a summary of the results. 



II. PROBLEMS IN PREVIOUS MODELS 

To construct our quasi-particle model, we first il- 
lustrate the problem hidden in early quasi-particle 
models. For simplicity, we take the scalar field as 
an example, similar problems also exist in the case 
of fields with spin. As is well known, one always 
uses the scalar field to characterize a system com- 
posed of Bose-type quasi-particlcs. Now, let us 
begin with the Lagrangian of quasi-particle with a 
temperature-dependent mass m = m{T) 

^ 1 1 

1=1 

where the effective mass m(T) describes the inter- 
action between real particles and N is the number 
of species of quasi-particles. Here, following Ivanov 
et.al [30|, we use this effective Lagrangian to de- 
scribe an interacting system composed of gluons, 
provided that we assume that there are only mas- 
sive and transverse quasi-gluons. 

Inserting the Fourier expansion of the field 0j {x) 



n= — oo p 



into this Lagrangian and according to the defini- 
tion of partition function in the path integral for- 



malism [3l|, |32| , one obtains 



N 



N'Yl 



;exp 



dT / d^xif 



Here the above formula expresses the partition 
function Z as a functional integral over (f>i of the 
exponential of the action in imaginary time and 
N' is an irrelevant overall normalization constant. 
Then, using the standard path integral procedure, 
we have the following partition function: 



InZ = VN 



l2nf 



1 



I3uj* - ln(I - e-'^'^ ) 



(2) 

where uj* = ^Jp^ + m? (T) is the dispersion rela- 
tion for the quasi-particle. From this expression 
we can see that the second term of the right- 
hand-side of Eq. (2) has the form of the ideal 
Bose gas formula for a quasi-particle except for the 
temperature-dependent dispersion relation coming 
from the effective mass rn{T); the first term of the 
right- hand- side of Eq. (2) is the zero point energy 
of the "vacuum" and is divergent. The difference 
between our case and that of standard statistical 
mechanics is that the infinity depends on the tem- 
perature T. In other words, the "vacuum" in our 
case is a thermal vacuum. In standard statistical 
mechanics we have m = const, and that divergence 
is independent of the temperature T, so we can 
throw this divergent part away in this case (this 
is because the vacuum zero-point energy and pres- 
sure cannot be measured experimentally and there- 
fore the zero-point energy should be subtracted). 
This operation will have no effect on computation 
of thermodynamical quantities. For example, by 
taking the derivative of the partition function, we 
can obtain the equation of state and the energy 
density: 



„ 1 dhiZ ^ d In,, 



')dp, 



1 d\nZ 
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-dp, 



where to = \Jp^^^t^rr? and d = N \s the degener- 
acy factor (d = 6 for SU(2); and d = 16 for SU(3) 
gluons). The above formula is formally model in- 
dependent because it comes completely from en- 
semble theory. We can see that the divergent part 
vanishes automatically after taking the derivative, 
because it is only an infinite constant. Now let 



us turn back to our current effective Lagrangian. 
When the function u!*{p) for the particle ("quasi- 
particle" in our case) excitation energy becomes 
temperature-dependent, the operation of taking 
the derivative with respect to (3 will no longer be 
valid. In particular, in the definition of energy den- 
sity, we find there is an extra term coming from the 
divergent part. If one wants to obtain physically 
meaningful results, he or she must treat this infin- 
ity carefully rather than throwing it away naively. 
A similar problem also emerges in the case of Fermi 
fields. For spin 1/2 field at zero chemical potential, 
from the Lagrangian of the quasi-fermion 



if = V^(z7^a^ - m(T))^, 



(3) 



one can do the derivation as before, and finally 
obtain: 



In Z = 2V 



(27r)3 



/3a;* + 21n(l + e-''" ) 



(4) 



where ut* = y^p'^ + m'^{T) is the dispersion re- 
lation for a Fermi-type quasi-particle. Similar 
to the case of Bose-type quasi-particle discussed 
above, there is also a temperature-dependent di- 
vergent term, which represents a temperature- 
dependent thermal vacuum energy. Indeed, just 
as was shown above, the two partition functions 
in both cases are ill-defined. If the infinity in the 
temperature-dependent thermal vacuum energy is 
not tackled properly, we will not be able to get 
physically meaningful results in the quasi-particle 
model. The aim of the rest of this paper is to 
treat this temperature-dependent "vacuum" en- 
ergy carefully and establish a thermodynamically 
consistent quasi-particle model. 



III. ELIMINATION OF THERMAL 
VACUUM DIVERGENCE 

In this section, we do our calculation for Dirac 
field at finite temperature and zero chemical po- 
tential and then turn to the case of scalar field. 
To eliminate the temperature-dependent infinity of 
the thermal vacuum energy, we introduce a classi- 
cal background field B (it is allowed to depend on 
the temperature) into the Lagrangian of the quasi- 
particle system. Thus, for a system composed of 
Fermi-type quasi-particles, we have the following 
Lagrangian 



if = ^{i-yi'df, - m{T))i^ + B. 



(5) 



Inspired by the approach of Walecka 133| , we shall 
use a series expansion method to separate the di- 
vergence term, and then choose a appropriate clas- 
sical background field B that satisfies the following 



condition 
B 



-Tr 



d^k 



\^ 



iv - 1)" 



(27r)4 Z^_^niji~ir 



(6) 



where 77 = — ^^ , as the counter term to remove the 
divergence and make the shift in the ground-state 
energy of the total system finite. 

After taking into account the effect of this clas- 
sical background field, by means of dimensional 
regularization we will have eliminated the diver- 
gence of vacuum zero point energy and obtain the 
finite, physically meaningful result of the shift in 
the ground-state energy 

Aeo= Eq — B — Eyac 



(47r)2 



m (T) In 
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mo 
13 



mo (mo - m{T)) 



mo (mo -m{T))^ 



-§(^0 



m{T)Y 



(7) 



This is what one expects in advance. Here we 
would like to stress that the series expansion 
method inspired by Walecka's approach for elimi- 
nating the temperature dependent infinity of vac- 
uum zero point energy given in Eq. (6) is the sim- 
plest way to achieve this goal. 

After completing the elimination of the infinity 
of thermal vacuum energy successfully in the Dirac 
field case, we turn our attention to the case of the 
scalar field. Similarly we start from the Lagrangian 
for Bose-type quasi-particle 



if = i9^(/.9'^0-i77i2(T)02 



B, 



(8) 



where 'rTi{T) is the effective mass of the quasi- 
particle. Using again the definition of measurable 
physical quantity and the series expansion, we can 
get the form of the classical background field 



B 



d^k 



It 



iv' - 1)" 

(27r)4"2^^^n(fc2_l)"^ 



(9) 



where rj^ 



%T) 



Owing to each boson's propa- 



gator having one more power of momentum in the 
denominator than that of the fermion, the classical 
background field B in a Bose-type quasi-particle 
system has only two terms. As is done previously 
for the Dirac field case, to calculate this integral, 
we generalize it to d dimensions and rotate to Eu- 
clidean space. We perform the integration and take 
the limit. In the end we get a finite shift in the 
ground-state energy 

1 



A£n = 



(47r)2 
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2 ^ ' mo 
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(10) 



Here we remark that the last term on the right- 
hand-side of Eq. (10) is a constant, and can be 
thrown away with no effect. Thus, with the help 
of the background field, we obtain a well-defined 
partition function of the quasi-particle system. 

The analysis of the divergence of vacuum zero 
point energy related to massless gauge bosons will 
be similar to the previous one, but the result is 
different compared to the case where the rest mass 
of the quasi-particle is nonzero. We repeat the 
above analysis and arrive at the form of the clas- 
sical background field for gauge bosons: 



B^N 



(27r)4 2 ^ 

^ ' n— 1 



^EM) 



-1 m^"(r) 
n(fc2-m2(r))' 



where TV is the degeneracy factor {N = 6 for SU(2); 
and A^ = 16 for SU(3)). Here we still adopt dimen- 
sional regularization to calculate E — E^ac and B. 
After some algebra we find: 



Aeq — Eq — Eyac — B 



N 1 



(47r)2 8 



-m\T). (11) 



Therefore, in the case of gauge bosons, there is 
only one term proportional to the fourth power of 
thermal mass contributing to the finite shift in the 
ground-state energy. 



IV. A SELF-CONSISTENT 
QUASI-PARTICLE MODEL 

We now use the quasi-bosons system as an exam- 
ple to illustrate the self-consistent statistical model 
for quasi-particles constructed above. It is well 
known that in equilibrium statistical field theory, 
once the partition function of the system is ob- 
tained, then all the thermodynamical variables can 
be determined. Therefore, the calculation of the 
system's partition function is very important. This 
is the reason why in this paper we take the parti- 
tion function as the starting point for studying the 
quasi-particle model. 

For Bose-type quasi-particle, the partition func- 
tion is 



lnZ= -y/3A£o-V" 



(2^ 



ln(l 



-Plo* 



), (12) 



where Aeg is the shift in the ground-state energy 
calculated in the last section for the scalar field. 
Then, according to the standard ensemble theory, 
we have 



1 dl-D-Z 

V = 

B dV 



-Aeo - T 



<fp 

(27r)3 



ln(l 
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dT 



Because the above formula is formally model inde- 
pendent, they automatically satisfy the fundamen- 
tal thermodynamical relation between the pressure 
ViT) and the energy density £{T) 

£{T)=T^^^-V{T). (13) 

For the same reason, we can generalize this formal- 
ism to the case of fields with spin. For the case of 
fermions, at zero chemical potential, the pressure 
and energy density are 



IdlnZ 

V = 

p dV 



1 d\nZ 

V dp 



-Aeo + 4T 



Aeo-T 



(fp 

l2^ 



5A£o 



ln(l- 



'), 
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dT 
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which follow from the partition function 



lnZ = -T//3Aeo+41/ 



d^p 

(27r)3 



ln(l 



-l3uo'\ 



(14) 

where Aeq is the effect of the classical background 
field in the quasi-fermions system. One can easily 
check that ^(T) and £{T) automatically satisfy 
the thermodynamical consistency relation (fT3|) . 

In Ref. 133| Walecka deals with the infinity of 
the vacuum zero point energy in the mean field ap- 
proximation at zero temperature. This is different 
from our case where the vacuum zero point energy 
is temperature dependent. Inspired by Walecka's 
approach, we have adopted the series expansion 
method to eliminate the temperature-dependent 
infinity of the vacuum energy at finite tempera- 
ture. From this procedure we have obtained a well- 
defined partition function and ensure that the ther- 
modynamical consistency relation is automatically 
satisfied. We also would like to point out that the 
above process of constructing quasi-particle model 
at finite temperature is general and can be directly 
generalized to the case of quasi-particle model at 
finite temperature and chemical potential. In this 
kind of system, the interaction between particles 
is related not only to the temperature but also 
to the chemical potential. In this case, we intro- 
duce a classical background field depending on the 
temperature and the chemical potential simultane- 
ously to eliminate the temperature and chemical 
potential dependent infinity of the vacuum zero- 
point energy. In this way, we can obtain a well- 
defined partition function. 



Here we would like to stress that a formal jus- 
tification or any direct relation to QCD is a real 
challenge for simp le quasi-particle models. The au- 
thors of Ref. [SJ] point out in their paper that the 
missing transport peak in the spectral function in- 
validates quasi-particle models for strongly inter- 
acting systems. On the other hand, the flexibility 
of quasiparticle models allows one to describe very 
well the thermodynamical variables of the quark- 
gluon plasma, as shown by various authors in var- 
ious examples J18l - l28l |. In this case, one needs to 
employ quasi-particle models to give predictions 
for more observables, as is done in Refs. [2^, [33 
(the authors of Ref. [331 adopted the equation of 
state obtained in the quasi-particle model to calcu- 
late the mass-radius relation of the quark star and 
find that the result is consistent with most recent 
astronomical observation). Through comparison of 
theoretical results with experimental data, one can 
constantly improve the quasi-particle models. We 
think that the best means to pursue this goal is 
to focus on this particular model, making connec- 
tions and comparisons with others whenever pos- 
sible and helpful. According to this thought, we 
use our model to fit the latest lattice results of S. 
Borsanyi et al for (2-1-1) flavor system [33 • Consid- 
ering the interacting plasma in thermodynamical 
equilibrium, we assume that it can be described 
as a system of massive quasi-particles, a picture 
arising asymptotically from the in-medium prop- 
erties of the constituents of the plasma. Just as in 
Ref. [21| , we write the effective mass of the quasi- 
particle as: 



mKT) 



n, 



where moi and 11^ are the rest mass and ther- 
mal mass of the quasi-particle, respectively. 11^ 
are given by the asymptotic values of the gauge- 
independent hard-thermal-loop self-energies: 



n„ 



■^qQ 



2ujgQ{mo + ujga) 



N' 



1, 



167V, 



T^G\T), 



n„ 



G^T) 



Nc + -Nf]T'G\T), 



{UN,-2Nf)\n{^) 



where Nc, Nf stands for the color factor and the 
number of flavors, respectively, and mg is the rest 
mass of the quark. T^/X is related to the QCD scale 
^QCD- The quantity G^ is to be considered as 



an effective coupling constant since it parametrizes 
all deviations of the exact spectral function from 
the "strict-quasi-particle" ansatz [2l| . Here wc use 
Tc = nOMeV as the transition temperature at 
vanishing chemical potential. Therefore, there is 
three free parameters: A, Tg and Bq (which is the 
integration constant in the pressure) in our quasi- 
particle model. By fitting the results of lattice 
QCD for the pressure V, the interaction measure 
{£ — SV) and speed of sound squared (Cf^)^ we 
can determine the following values of these param- 
eters: A = 6.6, Ts=-0.88Tc and B^^^ = SOMeV. 
As is shown in Figs. 1-3, our quasi-particle model 
can reproduce the lattice results well except the 
interaction measure {£ — 3V). According to the 
explanation of K.K. Szabo et al [2^, the reason 
for this can be easily understood. First, we use 
the pressure V in our fitting procedure, therefore 
a better agreement is expected for this quantity 
than for the interaction measure (S — SV), which is 
a prediction of the model. Second, the interaction 
measure can be obtained as a partial derivative of 
the pressure: 



S-3V 



dT\T^) 



A small difference in the pressure V can cause an 
evident deviation between the two results. In or- 
der to compare the results obtained in our model 
and those in earlier models, we also show the re- 
sults calculated according to K.K. Szabo model 
proposed in Ref. [2J]. From these figures, it can 
be seen that our model can fit the recent lattice 
data better than earlier models. Moreover, in their 
model, in order to fit the lattice result, the gluon 
degeneracy factor is fixed to be 16.4^0.2- So it 
is also a model parameter in earlier quasi-particle 
models. This is obviously unphysical. But in our 
model dg = 16, which is required by QCD. In other 
words, in our case there is no need to reserve this 
degree of freedom in the fitting procedure. 

Here we would like to discuss the main dif- 
ference between our model and previous quasi- 
particle models. In the previous quasi-particle 
models, in order to satisfy the requirement of thcr- 
modyiiamical consistency, those authors in Refs. 
|2(]| - [26l[28| have to introduce an additional medium 
contributions which they also call B(T). However, 
as discussed before, there exists a temperature- 
dependent infinity of thermal vacuum zero point 
energy and this makes the partition function of 
the previous models ill defined. Whereas in our 
model, by means of introducing a classical back- 
ground field, we have successfully eliminated the 
temperature-dependent infinity of thermal vacuum 
zero point energy. This makes the partition func- 
tion well defined and the thermodynamical consis- 
tency relation is automatically satisfied. Here it is 
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FIG. 1. Lattice data of (2+1) flavor QCD for the pres- 
sure P at T > Tc normalized by T" HI] (full boxes). 
The solid line corresponds to our quasi-particle model. 
The dashed line is the result obtained according to 
K.K. Szabo's model in Ref. [H 
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FIG. 3. Lattice data of (2+1) flavor QCD for the speed 
of sound squared (Cf^)^ in the QCD plasma [361] (full 
boxes). The solid line corresponds to our quasi-particle 
model. The dashed hne is the result obtained according 
to K.K. Szabo's model in Ref. [H. 
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FIG. 2. Lattice data of (2+1) flavor QCD for the in- 
teraction measure £ — SV at T > Tc normalized by 
T [3a] (full boxes). The solid line corresponds to our 
quasi-particle model. The dashed line is the result ob- 
tained according to K.K. Szabo's model in Ref. [24[ ]. 




FIG. 4. Additional medium contributions B normal- 
ized by r* in K.K. Szabo's model in Ref. [H (dashed 
line) and our classical background field B (solid line). 



also interesting to compare these two B{T) in the 
same graph for reference. The results are shown in 
Fig. 4. From Fig. 4 we can see that the overall 
trend of the curve of these two B agrees quali- 
tatively: as the temperature increases, they both 
decrease, although the rates are not the same. 



V. SUMMARY AND CONCLUSIONS 

To summarize, in the study of equilibrium sta- 
tistical field theory, it is quite important to estab- 
lish a well-defined partition function from which 
one can determine all the thcrmodynamical prop- 
erties of the system. The primary goal of this pa- 
per is to build a well defined partition function for 
a quasi-particle system whose effective mass de- 
pends on the temperature. We first adopt the se- 
ries expansion method inspired by J.D. Walecka's 
approach [33| to deal with the T-depcndcnt infinity 
of thermal vacuum energy by introducing a classi- 



cal background field into the effective Lagrangian. 
Through such a treatment we get a well-defined 
partition function. On the basis of this, follow- 
ing the standard ensemble theory we propose an 
improved quasi-particle model without any need 
of reformulation of statistical mechanics or ther- 
modynamic consistency relation. Then, as an ap- 
plication of our quasi-particle model we apply it 
to the case of (2+1) flavor QCD at zero chemical 
potential and finite temperature. It is found that 
compared to the earlier model, our quasi-particle 
model show a nice agreement with the latest lat- 
tice simulation results of S. Boranyi et al [36[. 
Moreover, in our model there is no need to reserve 
the giuon degeneracy factor as a model parameter 
in the fitting procedure, which is quite different 
with the previous models [2J]. Finally, We want 
to stress that our model framework is general and 
could be generalized to the case where the effective 
mass of the quasi-particle depends not only on the 
temperature but also on the chemical potential. 
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